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A LOCAL AND GLOBAL WELL-POSEDNESS RESULTS FOR 
THE GENERAL STRESS-ASSISTED DIFFUSION SYSTEMS 

MARTA LEWICKA AND PIOTR B. MUCHA 


Abstract. We prove the local and global in time existence of the classical solutions to 
two general classes of the stress-assisted diffusion systems. Our results are applicable in the 
context of the non-Euclidean elasticity and liquid crystal elastomers. 


1. Introduction and the main results. 

There are a number of phenomena where inhomogeneous and incompatible pre-strain is 
observed in 3-dimensional bodies. Growing leaves, gels subjected to differential swelling, 
electrodes in electrochemical cells, edges of torn plastic sheets are but a few examples 
da [Ml [351112]. It has also been recently suggested that such incompatible pre-strains may 
be exploited as means of actuation of micro-mechanical devices [MIETI. The mathematical 
foundations for these theories has lagged behind but has recently been the focus of much 
attention. While the static theory involving thin structures such as pre-strained plates and 
shells is now reasonably well understood [211 IHl EH El [27j , leading to the variationally re¬ 
duced models constrained to appropriate types of isometries mi El EH, and requiring bringing 
together the differential geometry of surfaces with the theory of elasticity appropriately mod¬ 
ified [2611301 I22j . the parallel evolutionary PDE model seems to not have been considered 

in this context. 


1.1. The model and the main results. In this paper, we are concerned with two systems 
of coupled PDEs in the description of stress-assisted diffusion. The hrst system: 


( 1 . 1 ) 


utt — {4>, Vu)^ = 0 

</.* = A(d^W(</.,Vu)). 


consists of a balance of linear momentum in the deformation field u : x M+ —)■ M^, and the 

diffusion law of the scalar field (/> : x M_|_ —)• M representing the inhomogeneity factor in the 

elastic energy density W. The field (p may be interpreted as the local swelling/shrinkage rate 
in morphogenesis at polymerization, or the localized conformation in liquid crystal elastomers. 

The second system is a quasi-static approximation of (jl.ll) , in which we neglect the material 
inertia uu, consistent with the assumption that the diffusion time scale is much larger than 
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the time scale of elastic wave propagation: 

— div 

( 1 . 2 ) 


,VM)j = 0 
cPs = A(^d^W{<P,Vu)y 


In both systems, the deformation u induces the deformation gradient, and the velocity and 
velocity gradients, respectively denoted as: 

F = Vu£ e rS, Q = V^i = Vv = Ft£ 


We will be concerned with the local in time well-posedness of the classical solutions to (|l.ll) . 
and the global well-posedness of (II.2p . subject to the (subset of) initial data: 

(1.3) m( 0, •)=tto, ut{0,-)=ui in R^, 


(1.4) = 4>o inR^, 
and the non-interpenetration ansatz: 

(1.5) detVrt>0 in R^. 
The main results of this paper are the following: 


Theorem 1.1. Let uq ~ id G ff'^(R^), ui G ff^(R^) and cpQ G i^^(R^). Assume that W is as 
in subseetion [23 Fix T > 0, and assume that the following quantities: 


( 1 . 6 ) 


|ui, Vuo - Ids, 


2 

H3 


/. 


-F ||uo - idllla + / W{(l)o,Vuo) dx 


are sufficiently small in eomparison with T, and with the eonstant 7 in |i. ?[ ). Then there 
exists a unique solution {u, (f>) of the problem 11.11) 11.31 - 11.51) . defined on the time interval 
[0,T], and such that: 


u - id G L°°(0, T; Uu G T°°(0, T; 

G L°°(0,r;ll'3(R3)) and fit ^ 


Theorem 1.2. Let fi^ G fl^(R^) and assume that W is as in subsection 11.31 Assume that 
||(/>o ||//2 is sufficiently small. Then there exists a unique global in time solution {u,fi) to 11.31) 
such that: 


u - id G T°°(R+; L®(R3)), V^n G L‘^{R+-, 

fi£ L°^{R+;H^{R^)) and Vfi £ L"^{R+-, (R^)). 


The proof of Theorem 11.11 relies on controlling the energy: 

[ + W{fi,Vu) dx, 

Jr3 2 

where the hyperbolic character of the hrst equation in (|l.ll) suggests to seek the a-priori 
bounds on higher norms of u and fi by the standard energy techniques. A detailed analysis 
reveals that the special structure of coupling in the stress-assisted diffusion system indeed 
allows for cancellation of those terms that otherwise prevent closing the bounds in each of 
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the two equations in (jl.ip alone. These terms are displayed in formulas (j2.6l) and (j2.8p in 
the proof of Lemma 12.21 Existence of solutions in Theorem 11.11 is then shown via Galerkin’s 
method, where we check that solutions to all appropriate e-approximations of the original 
system (jl.lj) still enjoy the same a-priori bounds in Theorem 12.31 This is carried out in section 
m while uniqueness of solutions is proved in section [H 

The proof of Theorem 11.21 given in section [U is based on the L^-approach as well. The 
system (II.2[) is of elliptic-parabolic type, thns there is no loss of regularity with respect to 
the initial data (in contrast to (jl.ljL . The analysis here is simpler than for (jl.lj) and we 
are able to show the global in time existence of small solutions. The toolbox we nse for 
the proofs of both resnlts is nniversal for hyperbolic-parabolic and elliptic-parabolic systems. 
Similar methods have been applied in [5l [121 ESI SSI ES] to study models of elasticity and 
their couplings with flows of complex fluids. A key element in these methods is the basic 
conservation law of energy and entropy type. 


1.2. The energy density W. We now introduce the assumptions on the inhomogeneous 
elastic energy density W in (II. ip . Namely, the nonnegative scalar field W : M+ x —>■ M_|_ 
is assumed to be in a neighborhood of (0, Ids) a-nd to satisfy, with some constant 7 > 0: 

W(0,ld) = 0, DIT(0,Id) = 0, and: 

(1 7) 

D^W(0,ld) : >j(l4>l^ + lsyniFl^) for all (0, F) G R x 


The two main examples of W that we have in mind, concern non-Euclidean elasticity and 
liquid crystal elastomers, where respectively: 


( 1 . 8 ) 


Wi(cf>,F)=Wo(FB(cf>)) + ^\c^l^, 

W 2 ( 0 ,F) = Wo(F(0)F) + i|<)-|2. 


are given in terms of the homogeneous energy density Wq : R^^^ R_|_ and the smooth 

tensor field F : R —)• R^^^. In both cases, we assume that i?((^) is symmetric and positive 
definite, and that B{0) = Id. Enrther, the principles of material frame invariance, material 
consistency, normalisation, and non-degeneracy impose the following conditions on Wq, valid 


for all F G R^''^ ^ ^ SO{3): 


(i) Wo{RF) = Wo{F). 

, . (ii) Wo{F) —)■ -boo as detF —)■ 0. 

^ (hi) ITo(Id) = 0. 

(iv) Wo{F) > c dist‘^{F,SO{3)). 

Examples of ITo satisfying the above conditions are: 


Wo,i{F) 

Wo, 2 (F) 


|(F^F)1/2 -Id|2 + llogdetFl? 


|(F^F )^/2 -Id |2 + 


1 


det F 


- I 


for det F > 0, 


where g > 1 and Wq,* is intended to be -boo if detF < 0 [3H]. Another case-study example, 
satisfying (i), (hi) but not (iv) is: Wo(F) = jF'^F — Idp. 
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We have the following observation, which we will prove in the Appendix: 

Proposition 1.3. For Wq which is in a neighborhood of SO{3) and B which is in 
a neighborhood of 0, assume fl.9\) and assume that B{Q) = Id. Then Wi and W 2 in Hl.^) 
satisfy 


1 . 3 . Background and relation to previous works. To put our results in a broader con¬ 
text, consider a general referential domain which is an open, smooth and simply connected 
subset of Let G : O —>■ be a given smooth Riemann metric on 17 and denote 

its unique positive dehnite symmetric square root by R = ^/G. The “incompatible elastic 
energy” of a deformation u of 17 is then given by: 

(1.10) E{u,n) = [ Wo{Vu{x)B{x)-^) dx Vu e 

Jn 

where the elastic energy density Wq is as in (|1.9I) . It has been proved in m that: 

inf E(u, 17) = 0 

uGVUl-2(n,K3) 

if and only if the Riemann curvature tensor of G vanishes identically in 17 and when (equiv¬ 
alently) the infimum above is achieved through a smooth isometric immersion u of G. 

It is worth mentioning that in the context of thin films when 17 = 17^ = 17 x (—§, with 
some 17 C there is a large body of literature relating the magnitude of curvatures of G 
to the scaling of inf E{-, 17^) in terms of the film’s thickness h, and subsequently deriving the 
residual 2-dimensional energies using the variational techniques. Firstly, in the Euclidean case 
of G = Ida, where the residual energies are driven by presence of applied forces /^ ~ /i“, three 
distinct limiting theories have been obtained [E] for iR(-,17'*) ~ h^ with j5 > 2 (equivalently 
a > 2). Namely: (3 € (2,4) corresponded to the linearized Kirchhoff model (nonlinear 
bending energy), /3 = 4 to the classical von-Karman model, and /3 > 4 to the linear elasticity. 
For /3 = 0 the membrane energy has been derived in the seminal papers 1S1125], while the 
case (3 = 2 was considered in M- Secondly, in |32j a higher order (inhnite) hierarchy of 
scalings and of the resulting elastic theories of shells, where the reference conhguration is a 
thin curved him, has been derived by an asymptotic calculus. 

Thirdly, in the context of the non-Euclidean energy (II.IUD . it has been shown in [7] that 
the scaling: inf ^E(-, 17^) ~ h^ only occurs when the metric G 2 X 2 on the mid-plate U can be 
isometrically immersed in with the regularity and when, at the same time, the three 
appropriate Riemann curvatures of G do not vanish identically; the relevant residual theory, 
obtained through T-convergence, yielded then a Kirchhoff-like residual energy. Further, in 
m the authors proved that the only outstanding nontrivial residual theory is a von Karman- 
like energy, valid when: inf \E{-, 17^) ~ This scale separation, contrary to [El [32], is due 
to the fact that while the magnitude of external forces is adjustable at will, it seems not to be 
the case for the interior mechanism of a given metric G which does not depend on h. In fact, 
it is the curvature tensor of G which induces the nontrivial stresses in the thin him and it 
has only six independent components, namely the six sectional curvatures created out of the 
three principal directions, which further fall into two categories: including or excluding the 
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thin direction variable. The simultaneous vanishing of curvatures in each of these categories 
correspond to the two scenarios at hand in terms of the scaling of the residual energy. 

Other types of the residual energies, pertaining to different contexts and scalings, have 
been studied and derived by the authors in [TOl [HI [HI Ell EH ESI EH EH [301 El] • 

Note that, at the formal level, the Euler-Lagrange equations of (ll.lOp are precisely the 
first equation in the system m- The dynamical viscoelasticity has been the subject of 
vast studies in the last decades (see for example IHEllIlElllElSI] and references therein), 
where various results on existence, asymptotics and stability have been obtained for a large 
class of models. For the coupled systems of stress-assisted diffusion of the type (11.11) . we 
found a substantial body of literature in the Applied Mechanics community [El [El [El EH 
EZ], deriving these equations from basic principles of continuum mechanics and irreversible 
thermodynamics. For example, the system derived in m is quite close to (jl.2l) from the view 
point of theory of PDFs; indeed the structure of nonlinearity in both systems is almost the 
same. However, derivation from the first principles aside, it seems that the analytical study of 
the Cauchy problem, particularly in long temporal ranges, has not been yet carried out. The 
closest investigation in this direction has been recently proposed in m, concerning existence 
of solutions for models of nonlinear thermoelasticity, and in |46j where the authors examine 
further models of thermoviscoelasticity from the viewpoint of mathematical well-posedness. 
We refer here to [El ES] as well. 


1.4. Notation. Throughout the paper we use the following notation. In (11.11) the operator 
div stands for the spacial divergence of an appropriate held. We use the convention that the 
divergence of a matrix held is taken row-wise. We use the matrix norm \F\ = (tr(F^F))^/^, 
which is induced by the inner product: {Fi : F 2 ) = iv{F'[F 2 ). 

The derivatives of W are denoted by DW, D^W etc, while their action on the appropriate 
variations (0, F) e R x R^^^ is denoted by: DW (0, Vtt) : (0, F), D‘^W (0, Vtt) : etc, 

often abbreviating to DW : {(l),F) and {D'^W) : (0, F)®^ when no confusion arises. 

The partial derivative of W with respect to its second argument is denoted hy dpW £ R^^^. 
The derivative in the direction of the variation F G R^^^ is then {{OpW) : F) G R. By 
(dpW) : (Fi (g) F 2 ■ ■ ■ ® Tfc_i) G R^^^ we denote the linear map acting on F G R^^^ as the 
kth derivative of W in the direction of Fi, F 2 ... Ffc_i, F. Hence, differentiating in F gives: 

: (Fi 8) ... Tfc) = ( {{d^pW) : (Fi ® ... Ffc-i)) : Ffc) G R. 

Finally, C,c> 0 stand for universal constants, independent of the variable quantities at hand. 


1.5. Acknowledgments. M.L. was partially supported by the NSF grant DMS-0846996 
and the NSF grant DMS-1406730. P.B.M. was partly supported by the NON grant No. 
2011/01/B/ST1/01197. 
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2. The crucial a priori estimate. 


Lemma 2.1. Every solution {u,(f)) to with regularity prescribed in Theorem M.ll sat¬ 

isfies: 

(2.1) 4 [ hut\^ + W{cpyu) dx+ [ |V(-A)-VtP dx = 0. 

dt Jr3 2 J 

Proof. Testing the first equation in (II. ip by ut and integrating by parts gives: 

4/ hut\^+W{fi,Vu) dx= f {ut,utt)+d^Wifi,Vu)fit + {dFW{ct),Vu):Vut) dx 
dt iR3 2 J 

= J {ut, div {dpWifi, Vu))) + {dpWifi, Vu) : Vut) + Vu)cj)t dx 

= [ d^w{fiyu)cj)tdx. 

Define = (—A)“^(/> and integrate the second equation in (jl.ip against fit- 

J dx = J = y 'fit^{dci>W{(p,Vu)) dx = - J ^d^W{(p,Vu)(j)t dx. 

Summing the above two equalities yields (12. ip and achieves the proof. ■ 

For every i,j, k £ {1, 2,3} we now define the correction terms: 

T^ijk —® VUxj. T ^Ux^^xis T ^Uxj^x^. ® 

+ {dldpd^W) : (V Uxi,Xj fixk H“ “t“ ^Uxj ,Xk4^Xi 

“t" V Xl-r . (b-r ■ Tt. ^■ Tl. H“ ^■ T ■ ) 


( 2 . 2 ) 


+ {d^d^pW) : Vur.. ® Vux, <8> Vu 

+ (50<92 IT) 

+ (5^5^ W/) 




{VUxi ® "^Uxjfixk + VUj;. (g) VUx^fixj + VUa: 

Ux^4>xjfix|^ T ^Uxjfixifixi; T 


Vu 


XkWxi 


) 


+ idt^)4>xi4>xj4>x^ 

+ iP(fy^){4’Xi,Xj4>Xk + fixpXkfixi + 4’xi,Xkfixj) ■ 

Lemma 2.2. Let (u,^) 6e a solution to /il.l\) . with regularity prescribed in Theorera M . 1[ For 
t > 0, define the two quantities: 


£it)= [ |ut|2 + |v3ut|2 + 2Ty((^,Vu) 
2K3 


, Vuj 


^Xi,Xj ,Xk^ ^ ^Xi,Xj ,Xk 


+ D^W{(j),Vu) : 

i,j,k=1..3 

Z{t) = ||'Wi|Ih 3(R3) + “ Id||^3(R3) + I|0||^3 (r3)- 

Then: 


\(g)2 


+ 2 J] nijkfix,,xpx, dx, 
i,j,k=1..3 


(2.3) 


< C(Z^ + z3/2). 
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Proof. 1. We differentiate the first equation in (jl.ip in a spacial direction Xi G {xi, 3 : 2 , X 3 }: 

Uxi,tt - div {dpW{(t),Vu) : = diY {dpd^W{(j),Vu)(l)xi) • 

We now differentiate the above twice more in the directions Xi^Xj G {xi,X 2 ,x^}: 

(2.4) Uxi^xj^xj^^tt div [dpW{cl),Vu) : Vuxi,xj,Xk) = div {dFd^W{(j),Vu)cj)xi,xj,Xk) +^i- 

The error term TZi above has the following form, where we suppress the distinction between 
different Xi,Xj,Xk, retaining hence only the structure of different terms: 

TZi = div(3(0|,Ty) : Vux ® Vuxx + 3{dpd^W) : {'Vuxxf'x + '^UxCpxx) 

(2.5) + 3{dpdlW)(fxf)xx + (d^W) : (Vuxf^ + 3{dld^W) : 

+ 3{d].dlW) : Vuxicfx? + {dpdlW)i(fxf^ . 


Above and in what follows, we also write {dpW) instead of dpW{4>,Vu), and {dpd^W) 
instead of d‘^d^W{4),'\/u), etc. Integrating (|2.4p by parts against Uxi^xj,xk,t we get: 


4 ^ ^ 

dt 


[ 

Jr3 2 


''Xi,Xj,Xk,t) d- (l)xi,Xj,Xk{{dFdcf)W) : '^Uxi,xj,xi,) 




( 2 . 6 ) 


/ 4^xi,xj,xk,t{i9p9(p^dd'j . ^Ux^ xj,xk} ds: 
JR3 


+ / </> Xi,Xj,Xl^ 

JM3 

T ~ I {dtdpdijfW^ . ^Uxi,xj,xp; ®^Ux^^xj,x}^ dx T f 'P'iUxi,xj,xi^,t dx. 

2 7r3 Jr3 


2. Differentiate now the second equation in (II.ip in Xi G {xi, 2:25 2 : 3 }: 

(fxa = A{{d^dpWi(f,Vu) : Vux,)+dlWi(f,Vu)(Px,). 

As before, differentiate twice more in Xi,Xj G {xi,X 2 ,X 3 }, to obtain: 

(2.7) 4>xi,xj,xk,t — 2^ {{d^dpW) : + id^W)(j)xi,xj,xk d-TZijk) > 

where TZ is given in ([^ . Testing ([221) against i-A)~^4>xi^xj,xk,t = i>xi,xj,xk,t, we get: 


( 2 . 8 ) 


/ \^'4>Xi,Xj,XkA dx — 

1 4>xi,xj,xk,t{{dFd(f,W) : ’S/ux^^xpxk) dx 

/r 3 

JR 3 

d 


-u 

{dt{dlW)){(l)xi,xpxu? dx + [ TZijkCpxi 
.3 iR 3 
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Note now that the hrst terms in the right hand side of both (12.61) and ()2.8p . namely the terms 
displayed in boxes, are the same. Consequently, subtracting (|2.8h from (|2.6h . we get: 


dx 


(2.9) 


1 d C 

“1“ / '^ijk4’xi,Xj,x^.,t dx 

Jr3 

— / 4^xi,xj,xi;{{dtOFd(f)W') . ^Uxi,xj,xi^} dx 

iR3 

+ “ / iptdpd^W) : '^Uxi^xj,xk ® '^Uxi^xj,xk dx 

+ “ Z' {dtid^W)){(f>xi xj,xk) dx + / T^iUxi,xj,xk,t dx. 
2 7r 3 ^ iR3 


3. We will now estimate terms in the right hand side of (j2.9p and prove that: 


( 2 . 10 ) 


and: 


f \(^x,,x,,xMdFd^W\\Vux,,x,,x,\dx+ f \dtd],d^W\\Vux,. 

Jr 3 Jr 3 




dx 


+ [ \dtidlW)\\4>x.,xj,xf dx<C (z3/2 + ||VV’i||^^3(M3)Z 


(2.11) I niUx,,x,,x,,t dx| < C (z3/2 + + C||VV’t|U3(R3) 




For the hrst term in (12.101) . we note that by the Sobolev embedding C'^d/2('j^3^ 
one easily gets: 

[ \(px,,XFxMdFd<f>W\\Vux,,x,,x,\ < \\{d].d^W) : Vut + {dFdlW)(^t\\L^\\V^cP\\L4^W\L2 

./R3 

< C (IIVUtllLoo + ||(/)t||Loo) ||V^())||/,2 II V^u||x,2 

< C(||Vut||^2 + ||A^t||^2) ||V3</,||^2||V"u||i2 
<c(z3/2 + ||V^t||^3(K3)Z). 

Similarly, the other two terms in (I2.10p are bounded by: 

Cdivuiiu- + Uth^) {W^Mh + Iiv^(^|li2), 

which implies the same estimate as before. 
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Regarding (j2.1ip . the first term in TZiUxi^xj,xk,t is bounded by: 

I |div(((9|.Ty) : Vux <8) Vuxx)\\'^^ut\ dx < c({\\Vut\\L°° + ||(/)t||i,oo)||V^n||i,oo||V^n||i 2 

Jr3 \ 

||V^'u||^4 “t" II7:^00 II^ II11^2 
< C (z3/2 + + ||VV’t||H3(R3)Z3/2) , 

becanse of the Sobolev embedding LP(M^) valid for any p G [2, 6]. Also: 

f |div((4fR) : (Vu,,)®3)||V\i| dx < c({\\Vut\\L^ + ||04||L-)||V2u||i6 
iR3 ^ 

+ llV^nlll^oo II V^'u||7;^2^ ||V^nj 117^2 

Other terms in TZi induce the same estimate as above. This establishes ()2.11l) . 

4. We now consider the last term in the right hand side of (12.9^ : 

( 2 . 12 ) f T^ijk4^Xi,Xj,Xk,t dx = ( ^ f T^ijk4^Xi,Xj,Xk f {T^ijk)t4’xi,Xj,X]^ dx. 

Jk 3 Vat 7^3 / 7 ]r3 

We now prove that: 

(2.13) I J^^{nijk)t(^x,,x,,x, dx| < C (z3/2 + + C||VlAi||773(E3) ^ ^ 

First, nsing the notational convention as in (12.5F TZ can be replaced by: 

TZ2 = 3{d^dpW) : Vux ® Vuxx + 2,{dldFd^W) : {Vuxx4’x + ^UxC^xx) 

(2.14) + {d^dlW) : {VuxT^ + KdldlW) : (Vu,)®Vx 

+ ^{dpdlW) : Vuxict^xf + {dlW){(t^xf + ^{dlW)(t>x(l>xx. 

The first term in (|2.14p can be estimated as before, using embedding and interpolation 
theorems: 


[ \{{d^dlW) : Vux Z Vu,^)J|V3(/.| dx 

iR3 
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V^llle + 
2 






Il 2 


while the third term in 7^2 is estimated by: 

7r3 

— C'^dlVirtlli.oo + ||(/>f||L' 

<C-(z5/2+Z2 + ||VlAdlH3(R3)Z 
Other terms in TZ 2 induce the same estimate as above. This establishes (I2.13p . 

5. Summing now (12.911 over all triples Xi,Xj,Xk, adding (12.111 . and taking into account 
(j2.12ll . (12.1011 . (12.1111 and (I2.13p . we obtain: 

P + {2\\VA\\l2 + llvVilli^) < o (z5/2 + z3/2) + C'||ViAdliu3(R3) {Z^ + Z) 

<d|ViAdlH3(R3) + c(z4 + z3/2). 

in view of Young’s inequality. Consequently, (12.311 follows and the proof is complete. ■ 

We now deduce the main a-priori estimate of this section: 

Theorem 2.3. Under the assumptions of Theorem 17.71 any solution on the time interval 
[o,r] to fop fll-Q satisfies: 

(2.15) supZ(7) <C'(T(0)+r2<fo(0) + ||no-id||2 2 ), 

where: To(0) = /jgs \ui\‘^ + 2W((/>o, Vuq) dx, and C is a universal constant. 

Proof. Assume that the quantities in (II.6p are sufficiently small. In particular, we require 
that Z{0) and that Z is sufficiently small on an interval [0,to]) where we also choose an 
appropriate to ^ T. Lemma 12.21 implies: E'{t) < CZ‘^/'^{t), which is equivalent to: 

(2.16) £{t)<c[ Z^/^{s) ds + £i0). 

Jo 

Further, by p2.ip it follows that: 

(2.17) sup||ut ||^2 < To(0). 

t 

Since: 

Vt < to ||tt(t) - id||^2 = 2 / f (n - id,ut) dx-k ||uo - id||^2 

Jo iR3 

< 

we easily obtain in view of (I2.17p : 


(2.18) 


S<t 


) ( sup u — i 

d|L2) 

s<t 

/ 


sup ||u(t) - id ||^2 < 4toSup ||ut(t )||^2 + 2||uo - id ||^2 

t<to 

<4t^To(0) + 2||uo-id||^2. 


(2.19) 
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[ 

Jrs 


Further, we observe that thanks to (|1.7jl . to Korn’s inequality and to Poincare’s inequality, 
there exist constants c, C > 0 so that: 
cZ{t) < 

[ \ut\^ + + 2W{ct),Vu) + V D^W{ct>,Vu) : 

+ [ |tt — idp dx < C'Z(t), 

ilR3 

as well as: 

^ dx| < c\\nL4^^nL^ < cz^/\t)z^/\t) = cz\t), 

i^j^k=1..3 

where we estimated each term in (|2.2I) by the Cauchy-Schwartz inequality and noted the 
appropriate Sobolev embedding. Consequently, we arrive at: 

(2.20) Vt < to £{t) + ||n - id ||^2 (t) > cZ{t) - CZ^{t) > cZ{t), 

provided that <C 1 is sufficiently small on the time interval we consider. In view of (j2.16jl . 
(j2.20p and (|2.19p . we now get: 

(2.21) Vt<to Z{t)<c(^J^ Z^/\s) ds + £{0)+tl£o{0) + \\uo-id\\l2y 

Calling Z = supjg[o,to] have: 

(2.22) Z< {CtoZ^/^'^Z + C{£{0)+tl£o{0) + \\uo-id\\l,), 

which combined with the requirement: CTZZ!'^ < 1 yields: 

(2.23) Z{to) <Z< C(^ (0) + tl£o{0) + ||no - id||2,). 

The above clearly implies the Theorem in view of the smallness of initial data in (jl.Op . ■ 

3. Proof of Theorem 11.11 Existence of solutions to fll.lh . 

In this section we construct approximate solutions to the Cauchy problem (11.11) (ll.3l - fm) . 
which satisfy the same a-priori bounds as in section [2l Given e > 0, consider the regularized 
problem: 


(3.1) 


utt — div^cIi?IT {(/): V^) j “ = 0 

(j)t = A(d^W{^,Vu)') 

with the same initial data as in doi-iiai). 

Lemma 3.1. Assume that all quantities in /il.6\) are sufficiently small. Then, there exists 
Te > 0 and a solution of h3.1\) U.A - \1.4^ on x [OjTe), such that: 

-id e L°°{0,T;H\R^)), ul^ £ L°°(0, T; ^^(M^)), 

(/>" G L°°(0,r;I73(M3)) and (fl e 
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Proof. 1. Since e > 0 is fixed, we drop the superscript ^ in order to lighten the notation in the 
next two steps. We proceed by the Galerkin method. Choose an orthonormal base 
in the space equipped with the scalar product: 

( 3 . 2 ) {'w,w)h4 = {w,w)l 2 + {'V'^w : V^w)l2. 

Similarly, let be an orthonormal basis in equipped with: 

( 3 . 3 ) {v, v)h3 = {v, v)l2 + ( V^u : V^v)l2 . 

Denote: = spanjrc^,..., and = spanju^,..., 

We now introduce the auxiliary scalar products: 


(3.4) 


{'w,w)w = {w,w)i 2 + (V^iu : V^w)i 2 \/w,'w € //'^(M^,M^), 
{v, v)v = (u, {—/S.)~^v)i 2 + (V^u : (—A)“^V^'i})i 2 Vu, v ^ H 


Clearly, these products are not equivalent to (13.2p . (13.3p . however their properties will allow for 
using the energy estimates of the proof of Lemma 12.21 to prove the regularity of approximate 
solutions {u^,4>^) which we define below. 

Let G X be the solution to: 


(3.5) 


« - dwdFW{(p^,Vu^) - eAu^,w^)^ = 0 
-Ad^Wicf^,Vu^),v^)^ = 0, 

ti^(0, •) = P^yiv(no), (u^)t(O, •) = Piuiv(tti), (O,-) ='^v^i4>o)- 


VZ : 1... A 

N! 


By P we denote here the orthogonal projections on appropriate subspaces. The classical 
theory of systems of ODEs guarantees existence of solutions to (j3.5p on some time interval 
[0,TAr). We now prove that these time intervals may be taken uniform for all sequences 

2. Since G and 4>^ G , (13.5p implies: 


(3.6) 


{utt - diYdpWicI)^,Vu^) - = 0, 

Note that the hrst equation in (13.61) is equivalent to: 

- diYdFW{^^,Vu^) - eAu^,uf 
+ {'^Z,XFXk,tt - dW{d],W{cj)^,Vu^) : - eAu^_^ 

^ (div(ai.a^W(</-^, + (7^f : V^uf), 


ij,k=1..3 

E ' 

ij,k=1..3 
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where by TZ^ we denote the error terms induced by the functions as in (j2.4p , (|2.5p . 

Likewise, the second equation in (j3.6p becomes: 


(cpf - Ad^W{4>^, Vu^), (-A)-Vf 




L2’ 




N 

Oui ^CC j ^cc ^ 


L2 


where we used the identity (j2.7p and the notation (j2.2p . with the superscript ^ indicating 
that they concern and 4>^. 

Let £[u ^be as in Lemma [22] with {u,(j)) replaced by and define: 

(3.7) £^[u^ + €{{-A)u’^ 


so that: 


[0^,u^](t)= [ |uf |V |V\fp+ 2kL((/>^,Vu^)+ e|Vu^P 

Jr3 


ij,k=1..3 

+ 2 ^ '^ijk^Z,xj,xk 

ij,k=1..3 




Xk I 


Following the proof of Lemmas 12.11 and 12.21 we find the counterpart of the inequality (12.3p : 

(3.8) £,[u^,(l)^]{t)<C f Z^/^[u^,(P^]{s) ds + £,{0), 

Jo 

where the constant C is independent from e, and where: 

Z[u^ ,4>^]{t) = ||llf^||^3(R3) + ||Vu^ “M||^3 (r3) + II</’^IIh3(R3)' 

Note that in order to obtain (13.8p we use only the equivalent formulations of (j3.6p above, 
hence indeed all the steps from the proof of Lemma 12.21 are valid with universal constants. 
Since the initial data in (13.5p consists of projections of the original data, their norms are 
uniformly controlled as well. 


3. We now consider the equivalence of with Z. Since for small Z one has: 

I n^Ax,,x, dx < CZ‘^l‘^[u^,<t>% 

we easily see that: 

(3.9) ,<p^] < CZ[u^ 

On the other hand, in view of (|3.7p : 

(3.10) £J(a^,4>^] > c,{Z - CeZ^/^) > c,Z[u^ 
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where by CejC'g we denote positive constants independent of N but depending on e. By (13. 8p 
we now arrive at; 

Z[u^ < Ce / ds + C'£fe(0). 

Jo 

Consequently, for to,e sufficiently small, we have: 

sup Z[u^ < Ce£e{0). 

t<tO,€ 

The above estimates, in particular (13.91) and (I3.1UD imply the uniform in N boundedness 
of the following quantities, on their common interval of existence [0, Tg]: 

- id G L°°(0, r,; uf G L°°(0, T,; 

yielding the weak-* convergence in as —)• oo (up to a subsequence), of the quantities: 

— id, , (p^, (p^ to the limiting quantities: — id, uf, 0^, (pf. Additionally, passing if 

necessary to a further subsequence and invoking a diagonal argument, we may also assure 
that: 

—)• Vix*^ and (p^ (p*^ point-wise in 

The Sobolev compact embedding: ^ C"((0,Te) x B{R)), valid on any 

ball B{R) C M^, justifies now that (p^ G C"((0,Te) x B{R)). Thus, in particular: 

(3.12) {d^W, dFW){(p^,Vu^) {d^W, dFW){^%Vu^) as iV ^ oo. 

It follows that {(p^,u‘^) is a distributional solution to (13.111 (|1.3I - [L^ . By (|3.11l) we obtain 
the desired regularity, completing the proof of Lemma 13.11 ■ 


Proof of Theorem 11.11 (existence part). Let cp^,u‘' be as in Lemma [3. II We first observe 
that a common interval of existence of {(p^,u^) can be taken as [0,T] with T prescribed by 
Theorem ll.il This follows through repeating the estimates in section[2l dealing with estimates 
of the first and the third order separately, and noting that the e-term appears exclusively in 
Te with a “good” sign. Consequently: 

(3.13) sup Zlu"^, (p^] < (17(^0, initial data), 

t<to 

and we see that indeed the solutions cp'^,u'^ can be extended over appropriate [0,T] with the 
quantities in (13.11) enjoying common bounds, independent of e. 

The same argument as in the last part of the proof of Lemma 13.11 implies now that the 
weak-* limit (up to a subsequence) of {(p'^,u'^) yield the desired regular solution {(p,u) to the 
original problem dni) ([L3]-[L1D. Condition (11.51) is automatically satisfied because of the 
smallness of initial data. ■ 
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4. Proof of THEOREM ll.lt Uniqueness of solutions to fll.ll) . 

Let {(j), u) and u) be two solutions to (II.ip with the same initial data. Define: 

{Scj)) = (j) — {6u) = u — u, 

and observe that: 

dFW{cl),Vu) - dFW{^,Vu) 

= dlW{^,Vu) : V{5u) + d^dFW{^,V^){6(j)) + D^dFWi^,Vu) : ((<5(/)), V((iu)) 
d^w{ci),Vu)-d^w{^,Vu) 

= dlW{^,Vu)i6cP) + d^dFWi^,Vu) : V{du) + D^d^W{^,Vu) : (((5</.), V(<5u))®", 


(4.1) 


where cj) and u are suitable linear combinations of (j), (p and u, u, given by the application 
of the Taylor formula. Subtracting equations (jl.ip for {4>,u) and {p,u) and using (14.ip . it 
follows that: 


{6u)u - div i^dlWip, Vu) : V{6u) + d^dFW{p, Vu 

= div (^D^dFW{P,Vu) : {{5(j)),V{6u))^'^ 

{5p)t - ^{dlW{p,Vu){6p) + d^dFW{p,Vu) : V(<5u)) 

= /\(D‘^d^W{P,Vu) : {{5cl>),V{5u)) 


We now test the first equation above by {6u)t-, while the second equation by (—A) to 

obtain: 

\ [ \6ut\‘^ + dlW{P,Vu) : (V((fu))®^ dx+ [ d^dFW{P,Vu) : {6P)V{Su)t dx 

^ JR3 J]r3 

= J D‘^dFW{p,Vu) : V(<5u)t) dx, 

\V{6p)t\^ + [dlWiP,Vu){5p)+d^dFWiP,Vu) : V(<5u)) dx 


(4.2) 


f D^d^WilVu) : {{6p),V{6u)f\6p)t dx. 


Consequently: 


hiU / dlW{(j),Vu){d(j))‘^ dx + f d^dFW{(j),Vu) : {6(j)t)'^{6u) dx 

^ dt J^3 ^ J]R3 

+ f \V{S<P)t\^ dx 


(4.3) 


f Vu) : {{6P), V(<5u)) dx 

iR3 


/. 


+ / dt 


Vtt)) ((5(/>)^ dx. 
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Adding (j4.2l) and (j4.3p . we arrive at: 

~ [ \{6u)t\‘^ + dlWilVu) : {Vi5u)f^ + dlW{^,Vumf 

+ 2d^dFW{(j),Vu) : {6^)V{5u) dx 
< C\\{6(j))t,(pt,ML^(R3)it) ■ srp 11 {6(p) , V (<5m) 11 ^2 (jrS) (t) , 


which implies that: 


(4.4) 


sup / + D^W{cP,Vu) : ({6(j)),V{6u))^^ dx < 

t Jr3 

Csup\\{6(p),V{5u)\\l2 [ \\{S(p)t,(pt,(pt\\L^{s) ds. 
t Jo 


As before, assumptions on W guarantee that the left hand side in (14.4p bounds from above 
the quantity: sup^ || V((^rt )|||2 • Since the integral quantity above is small for t <C 1, it 

follows by gaD that {5<p) and V{5u) are zero. ■ 


5. Proof of Theorem 11.21 The elliptic-parabolic problem (II.2|) . 

As in section [3l we first derive an a priori estimate for solutions of (|1.2p , whose existence 
will follow then via Galerkin’s method, in the same manner as for the system dni). 

Lemma 5.1. Assume that {4>,u) is a sufficiently smooth solution to which remains in 
a vicinity of (0, id) for all t >0, in the sense that: 

roo 

E[(j),Vu-ld] :=sup(||(/>||^ 2 (K 3 ) + ||Vu-/d||j:^ 2 )-bc IIV(/>, V^u||^2(r 3) dt < 1 
Then: 

sup(||Vu(t) -Id ||^2 + ll'u(t) -idll^e) + Vix - Id] < C'||0 o||h2(r3)- 

Proof. 1. Observe first the following elementary fact: 

(5.1) ||V^tt(t)||j^i(R3) < C'||V(/>(t)||jyi(R3). 

To see m, consider the first equation in (II.2p : 

(5.2) Vu) : Vu,, = -d^dpWicf, i = 1... 3. 

Condition (jl.7p and Korn’s inequality imply that the system (15.2p is elliptic, hence its solu¬ 
tions (normalised so that Vrt — Id E L®(M^)) obey: 

l|V^^lli,p(K3) < C'||V0||x,p(R3) p = 2,4. 

Differentiating (|5.2p with respect to x leads further to: 

liv^ulli^ < C{\\V^cf\\F2 + < C{\\V^cf\\F2 + ||V0||i4) < C{\\V^\\l2 + 

proving p5.ll) in view of the assumption in the Lemma. 
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We also observe the resulting control of pointwise smallness of (p and (Vtt — Id), by the 
Sobolev embedding: 


(5.3) 


sup \\<p, Vu — Idllioo < Csup \\4>, Vu — Id||j :^2 < <C 1. 


t>o 


t >0 


2. Testing the first equation in (II.2p by ut and the second one by ipt = (~^) ^(pt, we 
obtain the energy estimate, as in Lemma l2.II 

[ W{p,Vu)dx+ [ |V(-A)-Vt|^ da; = 0. 

To derive the second energy estimate we proceed slightly differently. Differentiating (jl.2p in 
a spatial direction Xi G {xi, X 2 , X 3 }, we get: 


(5.4) 


div (4iy(0, Vn) : Vn,, + dFd^WiP,Vu)P,,) = 0, 
= A{d^dFW{P,Vu) : Vu,,+dlW{(P,Vu)P,,). 


Now, testing the first equation above by Uxi, testing the second one by ^pxi = {—A) ^(pxi, 
and summing up the results, yields: 

1 d 


(5.5) 


2 dt 


\V{-A)-^Pxf dx+ f D^W{P,Vu) : {Px^'^UxT^ dx = 0. 

JR3 


Consequently, thanks to (HZD, the strict convexity of W implies: 

3 


2 dt J^3 


VVP + J [ (|V</.|2+^|(symVn,J|2) < 0 . 

2yR3 ^ 


Using Korn’s inequality and integrating in time we see that: 

(5.6) sup / (p'^dx + cf f (|V(/>|^ + |V^n|^) dxdf < C'||(/>o||^2fR3v 

t>o Jr3 Jo Jr3 


(5.7) 


3. We now differentiate (15.41) in a spatial direction Xj G {xi,X 2 ,X 3 }, getting: 

div {dj,W{(p,Vu) : Vuxi,xj + d^dFW{(p,Vu)(pxi,xj) = div7^l, 

(pxi,xj,t - A{d'^W{(p,Vu)(pxi,xj + d^dFW{(p,Vu) : Vuxi,xj) = A7^2, 

where the error terms TZi and TZi have the following structure (we suppress the distinction 
between different Xi,Xj): 

(5.8) 7^l,7^2 ~ D^Wi(P,Vu) : ((Vu,)®^ + {Vux)(Px + {(px?)- 

Integrating (15.71) by parts against Uxi^xj and {—A)~^<pxi,xj, respectively, it follows that: 


(5.9) 


[ \^P^x.,xf dx +c [ D^Wi(P,Vu) : {(Px,,x,,^Ux,,x,f^ dx 
dt Jr3 Jr3 

<C [ \V^u\^+ \V(p\^ dx, 
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because of (j5.8p and: 

I / {divTZi)uxi,xj dx|+| / (A7^2)(-A)■Vxi,XJ dx| 

iRS Jm.3 

< e\\Vuxi,xj, 4>xi,xj\\'i2(^3^ + C'||7^i,7^2||i2(R3)- 
Differentiating ()5.7p further, we obtain: 

div [dj,W{(j),Vu) : Vuxi,xj,xk + d^dFW{(l),Vu)(j)xi,xj,xk) = div7^3, 

(l>xi,xj,Xk,t - ^{dpV{(j),Vu)(pxi,xj,Xk + d^dFW{4>,Vu) : Vuxi,xj,Xk) = ^^4, 

where, as before: 


7^3,7^4 ~ D W (^, Vu) : ^V Uxx ^ ^^ '^xx4^x H" ^ 'Ux4^xx 4^x4^xx^ 

+ D^W{^, Vu) : ((Vn,)®3 + {Vux)®^<t^x + Vux{<t^xf + Vux{cl>x? + (0x)") • 


Testing by Uxi,xj,xk and 'tpxi,xj,xk, we find: 

[ \^‘^xi,xj,xu\ dx + c f D W{(j),\7u) : {4>xi,xj,xki^'^xi,xj,xk)^ dx 
inl dt JkS Jr3 

+ iv^ul^ + |V(/.|^ + |V20|4 + \v\f + |V(/.|® dx. 

Summing (|5.5p . (15.9p . (IS.lOp . integrating the result in time in the same manner as in (|5.6p . 
and recalling (|5.1I) . we obtain: 



POO P 

E[^, Vu - Id] < C / / |V\|^ + + iV^ul® + |V(/)|® 

Jo Jr3 

(5.11) + |V^u|^ + |V\|^ + dxdt + C\\J)o\\h2(m.3) 

POO P 

<c / \v<p\^ + \Vcj,f + |V20|4 dxdt + C||4 IIh2(m3)- 

Jo J 

We further have: 


POO P 

I / |V(/)|^ + |V(/.|® + |VVl^ dxdt 

0 7m3 

POO 

< csup {mwi^ + iivc^iii^ + / iiv</-i|2,, 3) dt < c(h2+h3) 

t>o Jo 


Conseouentlv. (15.111) becomes: ■=- < C{r? + + C'||(/)o||^ 2 - By the assumed smallness of 

E[(p,Vu — Id], we see that: 


(5.12) 


E<2CUo\\h2. 
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4. We now conclude the proof of the a-priori bound. Test (|5.2p by {u — id) to get: 
[ dlW{(l>yu) : (Vu - Id)^^dx 

<C f |<^||Vu-Id| + |M-id|(|Vu-Id| + |(^|)(|V(^| + |V\|) dx 

<C\ml. + {e + CE^/^)\\Vu-ld\\l, 

Indeed, by ()5.3p and (|5.1I) : 

[ \u- id| (|Vn - Id| + \cl)\) (| V())| + | V^n]) dx 

< C\\u - id||i6 ||Vn - Id||i 2 ||V(/), v20||^3 + C\\u - idU^e | 




< C\\Vu - Idlli^H^/^ + cml, + CE\\Vu - idlli. 

Thus, we obtain the bound on ||Vtt — Id||^ 2 , and subsequently on ||n — idU^e. 


A proof of Theorem II.21 Given {4>,u), consider the following problem which is the lin¬ 
earization of (11.21) at (0, id): 


(5.13) 


where: 

(5.14) 


div (^dlW{0,ld){Vu - Id) + d^dFW{0,ld)(p^ = divA, 

- A(5|lT(0,Id)(/. + 5,^5i.W(0,Id)(Vii-Id)) = AB, 

A = d^W(0, Id)(Vh - Id) + d^drWiO, Id)(^ - drWic^, Vu), 
B = d^W{^,Vu) - dlW{0,ld)^ + d^dFW{0,ld){Vu - Id). 


Let T be its solution operator, so that T[4>,u] = {(p,u). We will prove that T has a fixed 
point in the space X, where: 


A = {(0,11); (p e L^(R+-,H^{R^)), V(Vn - Id) e L°°{R+-, H\M.^)), 

V(j) e L'^{R+-,H^{R^))), V(Vn-Id) G L‘^{R+-, H^{R^)))}. 


Note first that the well-posedness of the system (|5.13l) follows by the Galerkin method in 
exactly the same manner as in section [3l under the regularity of the right hand side: 

A G L°°(M+; ff^(R^)) n L^(R+; ff^(R^)), B G L‘^{R+-, H^{R^)) 

Approximative spaces are constructed for cp G H^{R^) and for rt — id such that Vrt — Id G 
H^(R^). We leave this construction to the reader and note that it is simpler than the one for 
the system dni). As in the proof of Lemma 15.11 solutions to (I5.13P then satisfy: 

n _ rOQ r 

sup/ (p^ dx + y^ / / D‘^W{d,ld) ■. {(pxi■,VUa:^)^‘^ dxdt 

t>0 iR3 ^ Jo Jr3 

<C||VA,Vi?||i.(R^^^3(K3))+G||0o|li2, 
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which is obtained by testing with and (—A) Similarly, the second and third deriva¬ 

tives bounds eventually yield: 

sup II 011^2 (r3) + IIV*/*) — Id)||^2(K_^.^2(R3)) 

(5.15) 

< CIIVA, Vi?||^2(R_^.J^2(R3)) + C'II<('0|Ih2(R3))) 


while: 

sup ||Vn — Id||^2(R3) < C'll^llioo(R^;H2(K3))- 
Directly from (|5.14l) we observe that: 


ll^llL°°(R+;H'2(R3))f-,L2(R+;/f3(R3)) < C S[(/), Vu — H] 
l|•S||L2(R^;//3(R3)) < (7 H[(^, ViZ - Id] 


provided the quantity S is small. Then, by ()5.15l) : 

H[0, Vu - Id] < CH[0, Vu - Id]2 + Co||</.o||^2(k3). 

Based on the considerations from the part about the a priori bound we observe that: 

^[4>, Vn — Id] < 2 C'o||</>o||^2(']]j3)) 

provided that H is sufficiently small. Hence the operator T maps a ball B d X with a 
sufficiently small radius, into itself. Observe further that T is a contraction over B, whose 
fixed point yields the unique solution to the system (II.2p . Theorem 11.21 is proved. ■ 


6. Appendix: Proof of Proposition 11.31 


The first condition in ()1.7p is obvious. A direct calculation shows that: 

DWi{(l>,F) : ((^,F) = cl>^ + {DWo{FBm : ^FB'{^)) + {DAVoiBB^ : FB((/.)), 
which implies the second condition in (II.7p . Further: 

D2lTi(0,Id) : = \^\^ + D^^Woild) : (0H'(O))®2 

+ 2D'^Wo{ld) : (B'(0) ® F) + D'^W^iU) : F®^ 

= |(^|2 + D2lTo(Id) : {F + 


> |i^p + c sym F + (^B'(O) 


where we concluded from (|1.9p that DWo(Id) = 0 and that II^Ho(Id) is positive definite on 
symmetric matrices. We also note that: D^IT2(0,Id) = II^ITi(0,Id). To conclude the proof, 
it is hence enough to show that: 

(6.1) |(^p -|- |sym F + (^B'(O)p > c(|(^p -(- |sym Tp), 


for all (j) and F. Expanding the square in the left hand side, dividing by \4>\ and collecting 
terms, this is equivalent to: 

(1 - c + |H'(0)p) + (1 - c)|sym(iF)p 2(sym(iF) : B'{0)) > 0, 

9 9 










STRESS-ASSISTED DIFFUSION 


21 


which becomes: 

1 — c- - —sym(iF) H— ^ > 0. 

1 - C ^ (f) ' y/l - c 

The above inequality follows from: 1 — c — ■^^\B '> 0, which is true whenever c > 0 is 
sufficiently small. ■ 
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